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ABSTRACT 

C*~) ' We study motions of galaxies in galaxy clusters formed in the concordance ACDM cos- 

mology. We use high-resolution cosmological simulations that follow dynamics of dark 
matter and gas and include various physical processes critical for galaxy formation: 
gas cooling, heating and star formation. Analysing motions of galaxies and the prop- 
erties of intracluster gas in the sample of eight simulated clusters at z = 0, we study 
velocity dispersion profiles of the dark matter, gas, and galaxies. We measure the mean 
velocity of galaxy motions and gas sound speed as a function of radius and calculate 
qq ■ the average Mach number of galaxy motions. The simulations show that galaxies, on 

average, move supersonically with the average Mach number of « 1.4, approximately 
independent of the cluster-centric radius. The supersonic motions of galaxies may po- 
' tentially provide an important source of heating for the intracluster gas by driving 

weak shocks and via dynamical friction, although these heating processes appear to 
be inefficient in our simulations. We also find that galaxies move slightly faster than 
I 1 the dark matter particles. The magnitude of the velocity bias, b v « 1.1, is, however, 

smaller than the bias estimated for subhalos in dissipationless simulations. Interest- 
ingly, we find velocity bias in the tangential component of the velocity dispersion, 
but not in the radial component. Finally, we find significant random bulk motions of 
gas. The typical gas velocities are of order « 20 — 30% of the gas sound speed. These 
random motions provide about 10% of the total pressure support in our simulated 
k> clusters. The non-thermal pressure support, if neglected, will bias measurements of 

the total mass in the hydrostatic analyses of the X-ray cluster observations. 
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1 INTRODUCTION 

Clusters of galaxies are powerful probes of both the over- 
all cosmic evolution and evolution of galaxies in dense 
environments. Most of the cluster baryons are in the 
form of hot diffuse plasma and stars in galaxies, with 
a fraction of ~ 10 — 50% of stars in the intracluster 
medium. Interaction between galaxies and the hot intra- 
cluster gas is thought to be important in shaping proper- 
ties of both components. The compression of the interstellar 



medium (ISM) of an accreted galaxy 
can lead to a starformation burst ( 


and cluster tidal forces 
Dressier & Gundll983l 


Gavazzi et al.lll995l: iRose et alJl2001 


iKoopmann & Kennevl 


2004; Sakaie^ljj2002j]_or active galactic nuclei (AGN) ac- 



tivity (e.g., Miller & Owen 2003) and morphological trans- 
formation l|Moore et alJll996l : iKoopmann fc Kennevl 120041: 
iGnedinl l2003h . The subs equent removal of the ISM by 
ram pressure stripping jGunn fc Gotli Il972h and other 



z 



processes (e.g., iNulsenl Il982l) can suppress starforma- 
tion, dramaticall y affecting prope r ties of cluster galax ; 
population (e.g.. lOuilis et aJ J |200CA ISchulz fc Struck! 1200 
IVollmer et"al]|200lt I Vollmerl 1200317 ! Galaxies, in turn, may 
heat the intracluster medium (ICM) via supernova or 

AGN feedback and enrich it with heavy elements (e.g 
.. ■ ■ ii ■ ■ " ■■ ' n 



|KaiseiHl99ll: iMetzler ■ fc Evraidll994IValageas fc Silklll99' 

Churazov et al. 2002: Fabian et al. 2002; Scannapicco fc O 
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20041: 1 Voit fc Ponmanl2003URuszkowski et all2004lh Galax- 
ies can stir the surrounding gas as they move through 
the ICM driving turbulent flows and mixing the gas. 
The kinetic energy of galaxy motions can also be an 
impor tant source of the ICM heati ng (e.g., iMillerl Il98rj : 
lEl-Zant. Kim fc Kamionkowskilliool) . 

In the currently favoured Cold Dark Matter (CDM) 
model of structure formation, galaxies form in the extended 
dark matter halos as gas cools and condenses in the centre. 
After a galaxy is accreted by a cluster, the tidal forces can 



2 Faltenbacher et al. 



quickly strip the outer loosely bound regions of DM halos 
and can lead to a total disruption of the halo and, possi- 
bly, the stellar system. The dynamical evolution of galaxy- 
size dark matter halos in groups and clusters (often called 
the subhalos) has been the subject of many recent studies, 
which used a new generation of high-resolution dissipation- 



jGhiena et alJll998tlTormen et alJll998t iKlvpin et al||l999|: 


Colin et alJll99a lOkamoto fc Habelll99a IColm et al. 2000: 


Ghigna et alJ 2000 


iDe Lucia et alJ 120041: iKravtsov et al. 


2004 


: iDesai et alJ 


20041: iDiemand et al.1 hoo4 iGao et al. 


2004 


:lReed et alJl2004D. 



These studies find that abundance of subhalos is in 
reasonable agree ment with the observed abundance of 
cluster galaxies l|Moore et alJ Il99a IKravtsov et all 12004 
IDesai et all |2004 see, however, Diemand ct al. 2004 and 
IGao et all l2004h . although there are suggestive differences 
in the radial distribution (IDiemand et, all 12004 IGao et all 
20041 : lNagaj_fc_jCravtsoyl 120051) and circular velocity func- 



tions TDesa^t^lTl2004T) . Interestingly, velocity dispersion 
of subhalos is larger than that of dark ma tter by a factor up 
to ~ 1.4 in the inner regions of clusters (|C^lm_et_aL 2000; 
iGhigna et ai]l2000l : IDiemand et"al]|2004t IGao et al J 12001 . 
implying that dynamical estimates of cluster masses using 
galaxy velocity dispersion may be biased. However, apply- 
ing results of dissipationless simulations to the properties 
of cluster galaxies is subject to many caveats. Gas cooling 
during g alaxy formation increases central density of D M ha- 
los Ce.g- lBlumenthal et alJll98& iGnedin et al.ll2004l) . mak- 
ing them more resistant to tidal disruption. In the inner 
regions of galaxies stellar density is typically considerably 
higher than that of DM and the luminous component can 
survive the tidal forces eve n if the surround ing DM halo is 
completely disrupted (e.g.. lGao et alj|2004|) . It is therefore 
critical to study the properties of cluster galaxies forming in 
self-consistent cosmological simulations. 

Although a number of studies during the last decade 
used gasdynamics cosmological simulation s with cooling to 
study spatial distribution of galaxies fe.g..lKatz et alj |l992|: 
|Pearce^^ilJ ^99|: ^os hikawa et al] 1200 ll : IWeinberg et al] 
120041: IZheng et al.ll2004fi ~ only a few analyses directly ad- 
dresse d the motions of galaxies in cluster halos. iFrenk et al] 
( 1996) used simulations with cooling and starformation to 
study distribution and dynamics of galaxies in clusters. They 
found that the velocity dispersion of massive galaxies in the 
simulations is lower than that of the DM by ~ 20 — 30%, 
which they attributed to the or bital energy lo s s due to dy- 
namical friction. More recently, iBerlind et al] a2003fl anal- 
ysed mean motions of galaxies, identified as dense baryonic 
clumps, in halos of different masses in SPH simulations. 
These authors find that galaxies selected to have baryonic 
masses above a certain threshold move slower than DM in 
galaxy and group-size halos, the bias that disappears for 
more massive cluster-size halos. 

In this paper we present analysis of a sample of eight 
galaxy clusters formed in high-resolution cosmological simu- 
lations of the ACDM model. The simulations follow dynam- 
ics of DM and gas and include relevant cooling and heating 
processes and star formation. We use a combined sample 
of cluster galaxies, constructed using the friends-of-friends 
(FoF) algorithm, to study the average statistical properties 
of galaxy motions and compare them to the correspond- 



ing statistics calculated for the DM and intracluster gas. A 
complementary study of the radi al distribution of galaxies i n 
these simulations is presented in lNagai fc Kravtsovl i2005|) . 

The paper is organised as follows. In § |2] we introduce 
the simulation and describe the cluster sample which the 
subsequent analysis is based on. In § [3] we describe details 
of the galaxy finding algorithm and present analysis of the 
overall spatial and velocity distributions of galaxies. In § 2] 
we study the radial profiles of the velocity dispersion and 
velocity anisotropy for the DM, gas, and galaxies. The radial 
dependence of the average Mach number of galaxy motions 
in simulated clusters is presented in § 14.31 Finally, in §[^jwe 
summarise our main results and conclusions. 



2 SIMULATED CLUSTER SAMPLE 
2.1 Simulations 

In this study, we analyse high-resolution cosmological sim- 
ulations of eight group and cluster-size systems formed in 
the "concordance" flat ACDM model: f2 m = 1 — S1a = 0.3, 
Sl b = 0.021ft" 2 , h = 0.7 and a$ = 0.9, where the Hubble 
constant is defined as lOO/i km s _1 Mpc -1 , and as is the 
power spectrum normalisation on 8ft -1 Mpc scale. The simu- 
lations were done with the A daptive Refinement Tree ( ART) 
iV-body+gasdynamics code |K ravtsov 1 Klypin fc Khok hlovl 
ll997tlKravtsovll99a IKravtsov. Klypfn fc Hoffmanll2002h . a 
Eulerian code that uses the adaptive refinemen t in space and 
time and (non-adaptive) refinement in mass llKlvpin et al] 
120011) to reach the high dynamic range required to resolve 
cores of halos formed in self-consistent cosmological simula- 
tions. 

To set up initial conditions we first ran a low resolution 
simulation of an 80ft _1 Mpc box and selected eight clus- 
ters. The virial masses of clusters we selected range from 
r; 7x 10 13 to 3 x 10 14 /i _1 M . The perturbation modes 
in the lagrangian region corresponding to the sphere of five 
virial radii around each cluster at z = have then been 
re-sampled at the initial redshift, Zi = 49, retaining the 
previous large-scale waves intact but including addit ional 
small-scale waves, as described bv lKlypin et alj fcOOll) . 

High-resolution simulations were run using 128 3 uni- 
form grid and 8 levels of mesh refinement in the compu- 
tational box of 80/1 -1 Mpc, which corresponds to the dy- 
namic range of 128 x 2 8 = 32768 and peak formal resolution 
of 80/32, 768 ~ 2.44/i kpc, corresponding to the actual 
resolution of as 2 x 2.44 ~ 5ft -1 kpc. Only the region of 
~ 10/t -1 Mpc around the cluster was adaptively refined, 
the rest of the volume was followed on the uniform 128 3 
grid. The dark matter particle mass in the region around 
the cluster was 2.7 x 10 s hT 1 M Q , while other regions were 
simulated with lower mass resolution. 

As the zeroth-level fixed grid consisted of only 128 3 
cells, we started the simulation already pre-refined to the 
2nd level (I = 0, 1, 2) in the high-resolution lagrangian re- 
gions of clusters. This is done to ensure that the cell size is 
equal to the mean interparticle separation and all fluctua- 
tions present in the initial conditions are evolved properly. 
During the simulation, the refinements were allowed to the 
maximum I = 8 level and refinement criteria were based on 
the local mass of DM and gas in each cell. The logic is to 
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keep the mass per cell approximately constant so that the 
refinements are introduced to follow the collapse of matter 
in a quasi-lagrangian fashion. For the DM, we refine the cell 
if it contains more than two dark matter particles of the 
highest mass resolution specie. For gas, we allow the mesh 
refinement, if the cell contains gas mass larger than four 
times the DM particle mass scaled by the baryon fraction. 
In other words, the mesh is refined if the cell contains the 
DM mass larger than 5.42 x W 8 ^ 1 M = 2(1 - f b )m p or 
the gas mass larger than 1.81 x 10 8 /i _1 Mq = 4/(,m p , where 
m p = 3.16 x 10 8 /i _1 M and f b = Q b /^m = 0.1429. We 
analyse clusters at the present-day epoch as well as their 
progenitors at higher redshifts. 

Simulations included gasdynamics and several physi- 
cal processes critical to various aspects of galaxy formation: 
star formation, metal enrichment and thermal feedback due 
to the supernovae type II and type la, self-consistent ad- 
vection of metals, metallicity-dependent radiative cooling 
and UV heating due to cosmological ionising background 
llHaardt fc Madaull996l) . The cooling and heating rates take 
into account Compton heating/cooling of plasma, UV heat- 
ing, atomic and molecular cooling and are tabulated for the 
temperature range 10 2 < T < 10 9 K and a grid of metallic- 
ities, and UV inten sities using the Cloudy code (ver. 96b4, 
iFerland et al .1119981) . The Cloudy cooling and heating rates 
take into account metallicity of the gas, which is calculated 
self-consistently in the simulation, so that the local cooling 
rates depend on the local metallicity of the gas. 

Star formation in these simulations was done usin g 
the observationally-motivated recipe (e.g.. lKennicuttll998l) : 
P* ~ PgasA*i with t t — 4 x 10 9 yrs. Stars are allowed 
to form in regions with temperature T < 2 x 10 4 K and 
gas density n > 0.1 cm" 3 . No other criteria (like the col- 
lapse condition V ■ v < 0) were used. Algorithmically, star 
formation events are assumed to occur once every global 
time step Ato ~ 10 7 yrs, the v alue close to the observed 
timescales (e.g.. lHartmannll2002f) . Collisionless stellar parti- 
cles with mass m» = p»Ato are formed in every unrefined 
mesh cell that satisfies criteria for star formation during star 
formation events. The mass of stellar particles is restricted 
to be larger than min(5 x 10 7 h~ 1 M©,2/3 x m gas ), where 
m gas is gas mass in the star forming cell. This is done in 
order to keep the number of stellar particles computation- 
ally tractable, while avoiding sudden dramatic decrease of 
the local gas density. In the simulations analysed here, the 
masses of stellar particles formed by this algorithm range 
from « 10 5 to 7 x 10 s hr 1 Mq. 

Once formed, each stellar particle is treated as a single- 
age stellar population and its feedback on the surrounding 
gas is implemented accordingly. The feedback here is meant 
in a broad sense and includes injection of energy and heavy 
elements (metals) via stellar winds and supernovae and sec- 
ular mass loss. Specifically, in the simulations analysed here, 
we assumed th at stellar initia l mass function (IMF) is de- 
scribed by the iMiller fc Scalcl (Il979l) functional form with 
stellar masses in the range 0.1—100 Mq. All stars more mas- 
sive than M* > 8 M deposit 2 x 10 51 ergs of thermal energy 
in their parent cell 1 and fraction /z = min(0.2, 0.01M, — 
0.06) of their mass as metals, which crudely approximates 



1 No delay of cooling was introduced in these cells. 



id 


rsoo 


^500 


m dark 




ITLstar 


1 


606 


963 


1.111 


0.105 


0.091 


2 


661 


1047 


1.431 


0.137 


0.115 


3 


621 


988 


1.202 


0.120 


0.085 


4 


520 


830 


0.706 


0.074 


0.054 


5 


486 


771 


0.570 


0.043 


0.059 


6 


535 


851 


0.763 


0.078 


0.061 


7 


506 


806 


0.648 


0.061 


0.055 


8 


465 


743 


0.509 


0.044 


0.045 


mean 


550 


875 


0.867 


0.083 


0.070 



Table 1. Characteristic properties of the clusters. Column de- 
scription: (1) the identification number of each cluster; (2) the 
radius rsoo enclosing overdensity of 500 X p cr u in /i _1 kpc; (3) 
circular velocity at r 500 : v 500 = ^/GM 5Q0 /r 50 o in km/s; (4-6) 
the mass of dark matter, gas, and stars within rsoo in units of 
10 14 (i _1 Mq. The last row shows the mean values averaged over 
all eight clusters. 

the results of IWooslev fc Weaver! lll995l) . In addition, the 
stellar particles return a fraction of their mass and met- 
als to the surrounding gas at a secular rate mi oss — 
m, C (t - tbirth + Tp) - 1 with C = 0.05 and T = 5 Myr 
iJungwiert et al]l200lT) . The code also accounts for SNIa 
feedback assuming a rate that slowly increases with time 
and broadly peaks at the population age of 1 Gyr. We as- 
sume that a fraction of 1.5 x 10" 2 of mass in stars between 
3 and 8 M explodes as SNIa over the entire population 
history and each SNIa dumps 2 x 10 51 ergs of thermal en- 
ergy and ejects 1.3 Mq of metals into parent cell. For the 
assumed IMF, 75 SNII (instantly) and 11 SNIa (over several 
billion years) are produced by a 10 4 M stellar particle. 

2.2 Cluster sample 

Table lists properties of the eight clusters used in this 
study. The mass of each cluster is defined within the radius 
enclosing the cumulative density of 500p cr it, where p cr it is 
the critical density of the universe. This choice of overdensity 
is motivated by the fact that c lusters are on average more 
relax ed in their inner regions fevrard. Metzler fc Navarrol 
1996). Additionally, most of the X-ray radiation comes 
from region within rsoo so this radius is preferred in X- 
ray observations. In the following we use the radius rsoo, 
mass M500 = M(< rsoo), and circular velocity at rsoo, 
W500 = v&Msoo/rsoo, to normalise the physical cluster- 
centric radii, masses, and velocities. We will use the nor- 
malised quantities to obtain statistics averaged over all sim- 
ulated clusters. 



3 GALAXY SAMPLE 
3.1 Galaxy identification 

We use the frie nds-of-friends (FoF) algorithm (e.g., 
lEinasto et aflll984h to identify galaxies in the simulations. 
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Figure 1. The data points represent the properties of all the 
galaxies found in the eight simulated clusters at z=0 using the 
approach described in section 13. II The upper panel compares 
the stellar mass content within spheres of radius 50 h~ kpc and 
75 h~ 1 hpc. The data points which coincide with the dashed di- 
agonal belong to galaxies which contain all of their stellar mass 
within a radius of 50 h~ 1 kpc. The middle panel shows a compar- 
ison between the masses M;j obtained by FOF with 11 = 5kpc/h 
and the masses within a sphere of radius 50/i _1 kpc. Points lo- 
cated on the dashed diagonal represent 100% consistency of the 
masses obtained by the FoF approach and the masses gathered 
within a sphere of 50 h~ x kpc. The lower panel depicts the relation 
between galaxy mass and distance to the cluster centre. Where 
the distance of each galaxy is scaled by rsoo (see Table of the 
hosting cluster. 



The galaxies are defined as groups of stellar particles linked 
by the FoF with a certain linking length. We should note 
that the use of the FoF algorithm in dense environments 
can lead to misidentification of systems in two ways. First, 
distinct objects connected by an accidental bridge of small 
number of particles can be identified by the algorithm as 
a single object. Second, the algorithm can identify an un- 
bound group of particles which does not correspond to a 
real physical system. To minimise such problems, we adopt 
a two-step clustering analysis. 

In the first step, we identify groups with the FoF linking 
length of 11 = 5/i _1 kpc. This empirical value is similar to 
the typical sizes of galaxies and is well below the average 
distance between galaxies in a cluster. It is also close to 
the peak spatial resolution in the simulation. In the second 
step, we consider identified groups of more than 100 stellar 
particles and apply the FoF analysis with a smaller linking 
length of = 2 h~ 1 kpc only to particles belonging to these 
groups. For our analysis, we then consider only identified 
stellar particle groups (galaxies) with more than 50 particles 
linked with 11 — 2/i -1 kpc. 

Figure compares the stellar masses, M50 and M75, 
within the spheres of rso = 50 kpc and rrs = 75 /i -1 kpc 



radii around the identified stellar groups in all eight clusters. 
The upper panel of figure^shows that the dominant fraction 
of stellar mass resides within a radius of 50 h~ kpc. In a few 
cases the M75 substantially exceeds A/50. These rare cases 
correspond to close encounters between galaxies resulting in 
a significant overlap of their stars. This occurs most often 
when a small galaxy passes near the central cluster galaxy. 
The comparison of the FOF masses Ma and M50 shows that 
the FoF with 11 — 5/i _1 kpc correctly links all of the stellar 
mass in galaxies. The small deviation from the Mu = M50 
line at low masses is due to the contamination by the back- 
ground particles. The background of intracluster stars and 
stars associated with the central galaxies at small masses 
provides a small contribution to the mass within 50 h~ kpc 
for small sized objects. The contribution is evidently small 
however for our purposes and we did not attempt to correct 
for it. 

The lower panel of figure Q shows the cluster-centric 
distance of the individual galaxies (scaled by the 7-500 of the 
corresponding cluster) versus M50. There is a small defi- 
ciency of low mass (M50 < 5 x 10 9 /["'Mq) and high mass 
(A/50 > 2 x 10 11 h~ 1 M Q ) galaxies at small distances. The 
small mass galaxies are likely affected by the limited numer- 
ical resolution and tidal force of the cluster which leads to 
the tidal mass and possibly complete disruption. The min- 
imum particle number limit we imposed can then exclude 
such objects. At large masses, the deficiency is likely due to 
the efficient merging via dynamical friction. 



3.2 Galaxy number density profile 

In our analysis below we use the average number density 
profile of the compiled galaxy sample. The profile is shown 
in figure|5| For comparison we also show the average DM and 
gas density profiles. The profiles are calculated by averaging 
profiles of individual clusters with radii normalised to rsoo 
and densities normalised to the mean densities within rsoo- 
We fit the ga laxy number density profi le with the analytic 
NFW profile dNavarro et alJll99rl ll997h: 



"■galM = 



no 



(r/r s )(l + r/r s ) 2 



(1) 



The galaxy profile is well described by the NFW profile at 
r/rsoo > 0.5 but is flatter than the profile of DM at smaller 
radii. The difference may be caused by tidal mass loss and 
disruption experienced by galaxies in the core of the cluster. 
The difference is rather small, however. The overall galaxy 
distribution is rather similar to that of dark matter at the 
radii we reliably probe. The detailed analysis of the galaxy 
number density profiles in these simula tions and compar- 
isons with observations are presented in lNaeai fc Kravtsovl 
(2C 



3.3 Velocity distribution 

Before we proceed to the analysis of average galaxy motions, 
it is important to consider velocity distribution of cluster 
galaxies. If the velocity distribution follows the Maxwell- 
Boltzmann distribution both the velocity dispersion and the 
mean velocity are meaningful and are connected by the well- 
defined relation. Figure [3] shows the velocity distribution 
of the compiled galaxy sample and the Maxwell-Boltzmann 
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Figure 2. Average density profiles of galaxy number density 
(diamond symbols) and density of dark matter (solid line) and gas 
(dotted line). The error bars for galaxy number density indicate 
the 1<t confidence interval. The dashed line is a fit of the analytic 
NFW profile (eq. 0). The relative normalisation of the DM, gas, 
and galaxy profiles is arbitrary. 



distribution. The distribution was calculated using galaxy 
velocities normalised to the circular velocities, «5oo, of their 
host clusters. The Maxwell-Boltzmann distribution shown 
in the figure is given by 



f(v) = N 



3/2 



v 2 exp 



3v 2 



(2) 



where 030 is the three dimensional velocity dispersion. Here 
o"3D is calculated using normalised galaxy velocities, v/vsoo- 
Note that the dashed line in the figure is not a fit: it is the 
distribution given above calculated with the galaxy velocity 
dispersion measured for the galaxies in our compiled sample. 

The figure shows that the velocity distribution of the 
galaxy sample is well described by the Maxwell-Boltzmann 
distribution. One can therefore derive the mean velocity of 
the galaxies from their three dimensional velocity dispersion. 
We will use this relation below to compute the mean velocity 
of galaxy motions. 



4 GALAXY MOTIONS 

In this section, we analyse the motions of galaxies in the sim- 
ulated clusters. The analysis consists of three parts. First, 
we measure the mean radial velocity, v r , the three dimen- 
sional velocity dispersion, o~zt> , and its radial and tangential 
components, ay and a t , of dark matter, gas and galaxies di- 
rectly from the simulations. The errors for the galaxy mea- 
surements are computed using the a 2 as the basic values, for 
the 3D distribution as well as for the radial and tangential 
components. Then the 67% confidence interval is computed 
and by means of error regression transformed to the actually 
sh own error range o f the various a distributions (compare 
to lCoh'n et alj l200Cfi. Using the radial profiles of oy and at 
we can then compute the radial dependence of the veloc- 
ity anisotropy parameter, /3 (see eq. 0). We also measure 



Figure 3. The normalised velocity distribution of the compiled 
galaxy sample (diamonds) and the Maxwell-Boltzmann velocity 
distribution (dashed line for all galaxies within r < 5rsoo)- The 
error bars for galaxy distribution indicate the lo confidence inter- 
val. Note that the Maxwell-Boltzmann distribution is not a fit. It 
is calculated using the measured velocity dispersion of the galaxy 
sample. 



the sound speed c s of gas directly from the simulation. Sec- 
ond, we perform analyses using the equilibrium equations 
to model different components of galaxy clusters: the Jeans 
equation for galaxies and the hydrostatic equilibrium equa- 
tion for the intracluster gas. Using the above measurements 
as inputs, we use the Jeans equations to solve for the radial 
velocity dispersion profile of galaxies. Similarly, we apply 
the hydrostatic equilibrium equation to solve for the sound 
speed of gas. These results are then compared to the direct 
measurements from the simulations to assess the validity of 
the equilibrium equations in modelling kinematics of galax- 
ies and gas. Third, we study the average Mach number of 
galaxy motions as a function of cluster-centric radius. As we 
show below, we find that galaxies on average move superson- 
ically with the average Mach number of 1.3 — 1.4 throughout 
the cluster volume. 



4.1 Velocity dispersion, anisotropy, and sound 
speed profiles 

We measure the mean radial velocity, v r , and the three di- 
mensional velocity dispersion profile, azo(r) (and its radial 
and tangential components, ay(r) and at(r)) for dark mat- 
ter, gas and galaxies in radial bins centered on the minimum 
of cluster potential after subtracting the peculiar velocity of 
the cluster. The peculiar velocity of the cluster is calculated 
as the mass-weighted bulk velocity of dark matter enclosed 
within rsoo . We scale these measurements by the circular ve- 
locity W500 measured at rsoo of each cluster. After re-scaling, 
we compute the average profiles for the sample of eight sim- 
ulated clusters at z — 0. 

The upper four panels of figure[K]show the profiles mea- 
sured in the simulations. The net radial velocity is small 
throughout the cluster for both the dark matter and gas, 
except in the cluster core where there is a small net inward 
motion. Galaxies, on the other hand, show some net radial 
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velocity in some radial bins. However, the magnitude of the 
velocity, ~ 0.2«5oo, is small compared to the velocity disper- 
sion of the galaxies. Note that we find considerably smaller 
net radial velocity if we analyse each cluster at an epoch 
near 2 = 0, when it appears to be most relaxed. This indi- 
cates that the net radial velocity seen in the z — sample is 
likely due to incompletely erased galaxy groups within the 
clusters. 

We find that the three dimensional velocity dispersion 
of galaxies is biased with respect to that of dark matter. In 
other words, galaxies, on average, move faster than the aver- 
age speed of dark matter particles in clusters. The velocity 
bias, b v = o"gai/0"dm, is b v m 1.1 within the virial radius and 
disappears outside the virial radius of the cluster. Note that 
the velocity bias we find for the stellar systems is smaller 
than the velocity bias foun d for subhalos in dissipationless 
simulations, b v ~ 1.2 — 1.4 llCohr^e^alJI^OOOl: iGhigna et alJ 
l2000l: iDiemand et al.ll2004l : bao et alJl2004l) . Interestingly, 
we find that the velocity bias comes entirely from the tan- 
gential component of the velocity dispersion. As can be seen 
in figure|S] the radial velocity dispersion of galaxies and dark 
matter match well. 

Figure |K| also shows that gas within the virial ra- 
dius is not at rest but has non-zero velocity dispersion. 
We find that the 3D velocity dispersion of gas is ap- 
proximately constant at <7 gas ~ 0.4 — 0. 5^500 • In the in- 
ner regions of clusters gas moves with typical velocity of 
o"gas/\/3 ~ 0.2 — 0.3c s , but gas velocity dispersion be- 
comes ~ 0.5c s outside the virial radius reflecting the in- 
creasing strength of t he infall motions and lesser d egree 
of relaxation (see also |Naeai 1 Kravtsov fc Kosowskvl 2003; 
ISunvaev. Norman fc Brvanll2003r) . Note that these random 
velocities of gas contribute to the pressure support within 
clusters, in addition to the support from thermal pressure. 

The bottom-left panel of the figure |S] shows the radial 
profile of the velocity anisotropy parameter, 



75 kpc/h 



50kpc/h 



P(r) = 1 - 



(3) 



2a 2 (r) 

(see, e.g.. lBinnev fc Tremaine!ll987il . The figure shows that 
both the dark matter and gas tend to have a slight radial 
anisotropy. The anisotropy is nearly constant i n magnitude: 
/3 « 0.2 for dark matter (cp. lHoeft et alJl2004h and P « 0.4 
for gas. Note, however, the velocity anisotropy of gas drops 
within P ~ 0.2 in the centre. For galaxies, we find similar 
values of P on the outskirts of clusters. However, the value of 
p decreases gradually and reaches zero at r < 0.5rsoo. A fit 
to the act ual measurements in simulations using the fitting 
formula of ICoh'n et all feOOCh . 



P = Pn 



4r 



r 2 +4 



+ Po 



(4) 



gives best fit parameters of (p m , Po) = (0.76, —0.39). This fit 
is shown as a dashed line in figure |S] It describes the values 
measured in simulations quite well. Our best fit value of 0o is 
system atically smaller than the value found bv lColm et alJ 
(2000) in dissipationless simulations. This indicates that ve- 
locities of galaxies in our simulations have more enhanced 
tangential velocity dispersion compared to the subhalos in 
dissipationless simulations. 

The bottom-right panel in the figure|^|shows the average 
sound speed of intracluster gas. We compute this quantity in 



center of galaxy 






► 

V gal 


v gas 
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Figure 4. The sketch illustrates the measurements of the sound 
speed of gas in the vicinity of galaxies. We first identify the cen- 
tre of the galaxy (marked as the solid circle) and draw a cylinder 
with the diameter of 10/i -1 kpc aligned parallel to the direction of 
galaxy motion. We then measure the average sound speed of gas 
in a segment of the cylinder with the size of 50/i _1 kpc located 
75h~ x kpc ahead of the galaxy centre. To estimate the Mach num- 
ber for a galaxy, we use the relative velocity of the galaxy centre 
and the intracluster gas, v ga i — v gas , as a measure of the galaxy 
speed. 



two ways. First, we measure the average sound speed of gas 
in the vicinity of galaxies, as sketched in the figure |1J and 
average over the measurements of all galaxies in each radial 
bin. Alternatively, we measure the average sound speed of 
gas by simply computing the average of the sound speeds 
measured in individual mesh cells in each bin. The figure 
shows that the average sound speeds calculated in these two 
ways agree well. Note that the sound speed of gas is not con- 
stant but increases monotonically toward the cluster centre. 



4.2 The equilibrium analysis 

In this section, we apply the two equilibrium equations, the 
Jeans equation for galaxies and the hydrostatic equilibrium 
equation for intracluster gas to calculate velocity dispersion 
profiles under the assumptions that clusters are spherically 
symmetric and are in equilibrium. We then compare the 
results to the direct measurements of the velocity dispersion 
from the simulations to assess the validity of the equilibrium 
assumption. 



4-2.1 The Jeans equation 

We apply the Jeans equation to derive the average velocity 
dispersion of galaxies in galaxy clusters. For the relaxed sys- 
tems with no rotational support (y r = vg = 0), the Jeans 
equation is 



1 al 



n gal dr v—- ' r ' +2(3 f ~~^r~ (5) 

where P is the velocity anisotropy parameter in Eq. ©, 
n ga i denotes the number density of galaxies and $(r) = 
GMtot(r)/r is the gravitational potential of the cluster. In- 
tegration of the Jeans equation requires knowledge of (1) 
the total mass profile, M to t, (2) the number density profile 
of galaxies, n ga i and (3) the radial profile of the velocity 
anisotropy, P(r). In order to solve the equation for the sam- 
ple of clusters using the average profiles, it is convenient to 
work with the dimensionless variables. We, therefore, scale 
each variable using their values at r$oo and then construct 
an average profile for all eight clusters. For the number den- 
sity of galaxies, we use the best-fit NFW profile discussed in 
§ 13.21 (see also Fig.|5J. For the velocity anisotropy profile, we 
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Figure 5. The radial profiles of kinematic properties of dark matter (solid), gas (dotted) and galaxies (diamonds) , averaged over the 
sample of eight simulated clusters at z = 0. Errorbars show 1 a errors of the mean. Upper-four panels: the mean radial velocity v r 
(top-left), the three dimensional velocity dispersion a^jy (top-right), and its radial and tangential components, oy (middle- left) and at 
(middle-right) . The dashed lines in these panels are the results derived using the Jeans equation. Bottom-left: the anisotropy parameter 
fi. The dashed line is the best- fit j3 profile of galaxies based on the formula J1J. Bottom-right: the average sound speed of intracluster 
gas c a . The triangles are the measurements of sound speed of gas ahead of galaxies, as sketched in Fig.|l] The dotted line is an average 
of the sound speeds in individual mesh cells. The two dashed lines show the solutions of the hydrostatic equilibrium equation, where we 
assumed <r tur b = (thin dashed line) and 0"t ur b = CT 3D[gas] (thick dashed line). 



use the fit to the act ual measurements in simulations using 
the fitting formula of IColm et al.l (1200(1 see eq. @ above). 



Using the derived values of oy(r) and the best-fit /3(r) 
profile, we compute the tangential and total velocity dis- 
persion profiles, a t (r) and usd (r). These profiles are shown 
by the dashed lines in the corresponding panels of figure |S] 
All equilibrium solutions agree well with the values directly 
measured in simulations. 



4-2.2 Hydrostatic equation 

For a system in equilibrium, the hydrostatic equation 

1 dp gas _ d$ 
p g as dr dr 



(6) 



relates the gas density p gas and gas pressure p ga s, to the 
gravitational potential of the cluster, $(r). The gas sound 
speed is given by 



Pgas °turb 



(7) 
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Note that this expression takes into account the additional 
contribution to the pressure due to turbulent gas motions, 
where Oturb denotes the three dimensional velocity disper- 
sion of turbulent gas motions. 

Using the hydrodynamical variables and the total mass 
profile MtotO") measured directly from the simulations, we 
solve for the sound speed of gas using the Eq. @ and Q. 
To solve these equations, we perform similar scaling and 
interpolation discussed in the previous section on the mea- 
sured quantities. The dashed lines in the bottom-right panel 
of the figure [5] show the results of the calculation. The thin 
and thick dashed lines show the sound speed calculated as- 
suming er turb = and oturb = 0" 3D [gas] respectively, where 
<T3D[gas] is the velocity dispersion of the gas shown by the 
dotted line in the top-right panel. The figure shows that if 
random motions of gas are ignored, the sound speed is over- 
estimated by ~ 10%. This indicates that random motions 
contribute > 10% to the total pressure support in simu- 
lated clusters, which is in a greement with recent observa- 
tions JSchuecker et alJl2004|) . If we assume o" tur b = o"3D, the 
estimated sound speed of gas is in good agreement with the 
values measured directly from the gas density and thermal 
pressure. 



4.3 The average Mach number of galaxy motions 

To estimate the Mach number of galaxies, M = v/c 3 , we 
need to measure the average velocity of galaxy motions, v, 
and the average sound speed of the intracluster gas, c„. Be- 
low, we compute the Mach number of galaxies in two dif- 
ferent ways. First, we measure the average Mach number of 
galaxies directly from the simulation by measuring the veloc- 
ity of each galaxy and sound speed of gas immediately ahead 
of the galaxy in the direction of its motion, as sketched in 
figure^] Once the Mach number is measured for all galaxies, 
we compute the average Mach number in the spherical bins 
and average over all clusters in the sample. 

Alternatively, we can measure the average Mach num- 
ber of galaxies using the results of the Jeans and hydro- 
static equilibrium equations in the previous section. As we 
discussed in § 13.31 (see Fig. [3J , the velocity distribution of 
galaxies is well-described by the Maxwell-Boltzmann distri- 
bution. This means that the mean velocity of galaxies can be 
calculated using their three dimensional velocity dispersion 
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Figure 6. The Mach number of galaxies orbiting in the cluster 
environment. The diamonds show the direct measurements from 
the simulations, averaged over the sample of eight simulated clus- 
ters at z = 0, and the error bars indicate the la confidence in- 
terval. The dashed line is the result of the equilibrium equations, 
where we assumed that ot ur b = crgj) [gas] to account for the turbu- 
lent gas motions in the estimate of the sound speed of intracluster 
gas. The figure shows that galaxies move supersonically with the 
Mach number of 1.3 — 1.4 throughout the cluster environment. 



average Mach number of ~ 1.3 — 1.4 throughout the cluster 
volume. 

In the first approximation, the slightly supersonic mo- 
tions of galaxies can be explained by the following simple 
argument. The mean velocity of galaxies is given by eq. (J8J, 



.1/2 



while the sound speed is v a = (7-P gas /p) = (7<r gas /3) 
where <7g as is the velocity dispersion of the microscopic ther- 
mal motions of gas particles. If the three-dimensional ve- 
locity dispersions of gas particles and galaxies are approxi- 
mately equal, as is expected in equilibrium, we have: 



(M) = v/c a « 0.92^/3/7 « 1.23, 



(9) 



for 7 = 5/3. As was shown in the previous section, turbulent 
pressure contributes > 10% to the total pressure and the 
sound speed is actually « 10% lower than what it would be 
without the turbulence. This can explain the average Mach 
number that we measure (M) « 1.23/0.9 ~ 1.37. 



0.92cr 3 D. 



(8) 



We thus compute v from 030 obtained by solving the Jeans 
equation. Similarly, we use the sound speed of gas, c s , ob- 
tained from the hydrostatic equation. Note that estimating 
the sound speed of gas requires the knowledge of turbulent 
gas motions uturb- Since characterising the nature of the 
turbulent gas motions is beyond the scope of the paper, we 
simply assume a tur b = C3D[gas], which has been shown to 
result in good agreement between the calculated and mea- 
sured values of the sound speed in Fig. [£] 

Figure E| shows the average Mach number of galaxies in 
the simulated clusters. The diamonds are the direct mea- 
surement from the simulation, and the dashed line is calcu- 
lated using solutions of the equilibrium equations. The figure 
shows that on average galaxies move supersonically with the 



5 DISCUSSION AND CONCLUSIONS 

We presented analysis of the galaxy motions in clusters using 
high-resolution cosmological simulations of the concordance 
fiat ACDM model. The simulations follow dynamics of dark 
matter and gas and include various physical processes crit- 
ical for galaxy formation: gas cooling, heating and star for- 
mation. These simulations, therefore, follow the formation 
of galaxies and their evolution in the dense cluster environ- 
ment in a realistic cosmological context. Analysing motions 
of galaxies and the properties of intracluster gas in the sam- 
ple of eight simulated clusters at z = 0, we study velocity 
dispersion profiles of the dark matter, gas, and galaxies. We 
measure the mean velocity of galaxy motions and gas sound 
speed as a function of radius and calculate the average Mach 
number of galaxy motions. 
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Our simulations show that galaxies, on average, move 
supersonically throughout the cluster volume. The average 
Mach number is (M) « 1.4, approximately independent of 
cluster-centric radius. The value of (M) can be attributed 
to the difference between the three-dimensional velocity of 
galaxies and one-dimensional sound speed and to the exis- 
tence of the turbulent motions of gas (see eq. |5J|. The ther- 
mal pressure and, hence, the sound speed are smaller than 
would be required if all of the gas pressure support was due 
to thermal pressure. Also, gas and galaxies have somewhat 
different radial density and velocity anisotropy profiles. The 
motions of galaxies and gas particles are governed by the 
same potential, but non-zero velocity anisotropy and differ- 
ences in the radial distribution can lead to different velocity 
dispersion and sound speed profiles. 

We also find that galaxies move slightly faster than the 
dark matter particles in clusters, although the magnitude 
of the velocity bias, b v « 1.1, is considerably smaller than 
the bias estimated for subhalos in dissipationless simulations 
JColm et al]l2000l:lGriigna et aUESoti iDiemand et al.ll2004t 
iGao et al.ll2004h Interestingly, we find velocity bias in the 
tangential component of the velocity dispersion, but not in 
the radial component. Despite the small sample size used 
in this analysis, we find that these results are robust and 
statistically significant. Namely, the results do not change 
if we analyse clusters by removing two of the least relaxed 
clusters at z = or using clusters at the epoch when they 
appear most relaxed. Nevertheless, it would be good to verify 
these results with higher resolution simulations and better 
statistics. 

Our simulations show that the difference between the 
dark matter and galaxy velocity dispersions is significant 
only in the inner regions of clusters (r/r^oo < 0.7). This is 
also where we find a corresponding difference in the veloc- 
ity anisotropies. With these results in mind, we conjecture 
that the bias is likely caused by the preferential disruption 
of objects on highly radial orbits. In other words, galaxies 
on radial orbits are prone to more efficient tidal disruption 
than those on circular orbits with the same orbital energy, 
since galaxies on radial orbits pass closer to the dense cen- 
tral regions during their peri-centric passages. If the effect 
is significant, we expect the fraction of galaxies with the 
small (large) tangential (radial) velocity component to de- 
crease in the inner regions of clusters. The distributions of v r 
{vt), therefore, become progressively skewed towards smaller 
(larger) values with decreasing radius. 

Our simulations show that this is indeed the case. 
Specifically, we find that the distributions of v r and vt be- 
come increasingly non-gaussian, and Vt/v r ratio increases 
monotonically at the smaller cluster-centric distances within 
v < rsoo. Furthermore, if galaxies accrete onto clusters with 
the positive orbital velocity anisotropy, j3{r) > (i.e., pref- 
erence for the radial orbits), similar to dark matter parti- 
cles, the tidal disruption would also drive the system toward 
more isotropic orbits. We indeed find that galaxy motions 
in the central regions of the clusters are nearly isotropic 
(/3 « — 0.1), while the DM particles have the radial 
anisotropy of /3 ~ 0.2 — 0.3 (see Figure [KJ. 

We also find considerable random bulk motions of gas. 
The 3D velocity dispersion is approximately constant as a 
function of radius: <x gas ~ 0.4 — 0.5wsoo- In terms of the 
sound speed, the gas moves with the typical velocity of 



a eas /V3 ~ 0.2 - 0.3c a in the inner regions of clusters (see 
also iNaeai et alJl2003l : ISunvaev et al J 12003). Outside rsoo, 
the typical velocities are ~ 0.5c 3 , which reflects the increas- 
ing strength of the infall motions and lesser degree of relax- 
ation. The random motions of gas contribute to the pressure 
support of galaxy clusters, in addition to the support from 
thermal pressure. We show that in our simulations, random 
motions contribute ~ 10% of the total pressure support. 
Recent X-ray observations of the Coma cluster also show 
evidence of random m otions of gas of a similar magnitude 
JSchuecker et alJl2004l) . The non-thermal pressure support, 
if neglected, will bias measurement of the total mass in the 
hydrostatic analyses of the X-ray cluster observations by 
~ 10 - 20%. 

The supersonic motions of cluster galaxies may be an 
important source of heating of the intracluster gas. Su- 
personically moving galaxies and group s can drive weak 
bow s hocks and shock-heat th e ICM (e.g. iMarkevitch et all 
120041: iFinoeuenov et al.1 12004) or deposit energy via vis- 
cous dissipation of soundwaves iRuszkowski et al]l2004l) . It 
is likely, however, that a more efficient heating mechanism 
is the transfer of the orbital energy of galaxy motions to 
the internal energy of gas via dynamical friction. Dynam- 
ical friction in the gaseous me dium is efficient only if the 
perturber moves superson i cally {Rude rman fc Spie^e]||l97ll: 
iRephaeli fc SalpeteJlT98ci: lJust et alJll990t IOstrikerlll999l) . 
Our results therefore co nfirm that this is a potentially vi- 
able heating mechanism. lEl-Zant et alj (I2004T) have recently 
argued that the dynamical friction heating is sufficiently 
efficient to prevent formation of the cooling flows in clus- 
ter cores. Their estimates show that in the Perseus clus- 
ter galaxies do move supersonically with the average Mach 
number of M ~ 1.5 — 2. For reasonable choices of the galaxy 
mass-to- light ratios the rate of energy loss to dynamical fric- 
tion is sufficient to offset a cooling of gas in this cluster. 
lEl-Zant et al.l i2004l) also emphasise that this mechanism is 
self-regulating: as gas is heated, the galaxy motions become 
subsonic and heating becomes inefficient. 

In principle, dynamical friction heating should be mod- 
elled self-consistently in cosmological simulations. Our sim- 
ulations, however, do suffer from the well-known "overcool- 
ing problem:" the fraction of baryons in stars and cold gas 
is at least a factor of two higher than observed for the sys- 
tems of the mass range we consider. For example, the typical 
star formation rate in the cool cores in the central galaxy 
in our clusters is > 1000 MQyr" 1 , which is much larger 
than is allowed b y the observed values of the mass accre- 
tion rates (see e .g.,|Kaa stra et all2004l : |peterson et all200ll . 
120031: iTamura et alJl200lD . This indicates that the dynami- 
cal friction heating is not efficient in our simulations. This 
could be due to their limited spatial resolution (a few kpc). 
It is possible, that to resolve the dynamical friction wakes 
properly a sub-kiloparsec resolution is needed. However, it 
is also possible that the wake formation is prevented by the 
random motions of the gas that we find in our simulations. 
These motions are much stronger than the gravitational pull 
of any individual galaxy and it is likely that galaxies simply 
cannot form wakes in such highly chaotic velocity field. 

Given that dynamical friction can provide an attractive 
source of the ICM heating, it will be important to pursue this 
subject further. The progress can come both from the higher 
resolution cosmological simulations and from controlled gas- 
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dynamics experiments of a gravitating body moving in a gas 
flow. It would be interesting, for example, to test whether 
dynamical friction is equally efficient for a perturber moving 
in laminar and strongly turbulent flows. 
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